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In conventional superconductors, magnetic impurities form an impurity band due to quantum 
interference of the impurity bound states, leading to suppression of the superconducting tran- 
sition temperature. Such quantum interference effects can also be expected in d-wave super- 
conductors. Here, we use scanning tunneling microscopy to investigate the effect of multiple 
non-magnetic impurities on the local electronic structure of the high-temperature superconductor 
Bi2Sr2Ca(Cui_i;Zni:)208+i . We find several fingerprints of quantum interference of the impurity 
bound states including: (i) a two-dimensional modulation of local density-of-states with a period 
of approximately 5.4 A along the a- and 6-axes, which is indicative of the d-wave superconducting 
nature of the cuprates; (ii) abrupt spatial variations of the impurity bound state energy; (iii)an 
appearance of positive energy states; (iv) a split of the impurity bound state. All of these findings 
provide important insight into how the impurity band in d-wave superconductors is formed. 



I. INTRODUCTION 

Impurity atoms play an important role in solid state 
physics, manifested for example by the significant 
doping dependence of the electrical conductivity in 
semiconductors. The effect of impurities is also crucial 
in superconductors. In conventional superconductors, 
magnetic impurities act as pair breakers and thus give 
rise to a suppression of the superconducting transition 
temperature (Tc), which therefore depends on the 
impurity concentrationiiS,. When a number of magnetic 
impurities are present, the wave functions of electrons 
bound at each impurity potential interfere with each 
other, leading to the formation of an impurity band with 
a width that is intimately connected to the reduction of 



states associated with each impurity site interfere with 
each other in the same manner as for conventional 
superconductoraii"— . Some theoretical studies have 
predicted that this interference leads to a splitting and 
shift in energy of the bound state in the LDOS, as well 
as to formation of the impurity band, that is closely 
related to the impurity effect on the bulk propertiesi^Ti^i. 
However, despite a great deal of theoretical effort, few 
experimental studies have investigated the effect of 
multiple impurities on the LDOS. 



In d-wave superconductors such as the cuprates, both 
magnetic and non-magnetic impurities suppress the 
bulk T(^i^ and modify the local electronic structure^rJ^. 
The local effect of such impurities has been studied 
using scanning tunneling spectroscopy (STS), reveal- 
ing that impurity bound states are created within 
the superconducting gap, which is (not) suppressed 
around non-magnetic (magnetic) impurity, in the local 
density-of-states (LDOS). For a single non-magnetic 
impurity such as zinc (Zn)^iiSi^ii^, the bound state 
lies close to the Ey as a near-zero-energy peak (NZEP) 
and the superconducting gap value is suppressed at Zn 
site. In real space, the bound state extends into the 
region surrounding the impurity site and adopts d-wave 
symmetry, referred to as a " cross-shaped pattern" . 
This pattern is characterized by maximum intensity at 
the impurity site and secondary maxima at the four 
second-nearest-neighbor sites^iii^ii^ii^. 

When a large number of non-magnetic impurities 
are present in a d-wave superconductor, it is theoret- 
ically anticipated that the d-wave symmetric bound 



II. EXPERIMENT 



Here, we use scanning tunneling microscopy (STM) 
to study the effect of multiple impurities on the LDOS 
of single crystal Bi2Sr2Ca(Cui_2:Zn2;)208+5 (Zn-Bi2212, 
with Tc — 87K for the as-grown sample). The single 
crystal used in this study was grown by the floating zone 
method. We chose a relatively high Zn concentration of 
X ^ 0.5% (the nominal value of x used in the crystal 
growth was 2%), which was determined by inductively 
coupled plasma optical emission spectroscopy performed 
on the as-grown sample. The STS measurements were 
performed in a helium gas environment at 4.2 K on a 
surface that was atomically clean and flat [see the inset 
of Fig. IHa)]; the surface was prepared by cleaving the 
sample in-situ at 4.2 K in pure helium gas. We mea- 
sured the I-V characteristics and obtained the tunneling 
spectra G'(r, Vb) = dl/dV by numerical differentiation. 
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FIG. 1. (Color online) Typical topographic STM image (field 
of view 3.5 x 3.5 nm^) of a BiO plane of Zn-doped Bi2212; the 
Bi atoms are visible as bright spots. The image was taken at 
a sample bias voltage of Vb = +300 mV and a tunneling cur- 
rent of Jsct = 400 pA. (b) Conductance G(r, Vb) map at a bias 
voltage of Vb = -3 mV for the same field of view as in (a). The 
bright central spot in (b) coincides with the Zn site indicated 
by the yellow circle in (a), (c) Schematic representation of 
the Cu02 plane. The intensity maxima (bright spots) in the 
conductance map of (b) correspond to the Zn atom, second 
nearest-neighbor (2nd NN) and third nearest-neighbor (3rd 
NN) sites, as indicated by the red (bold) circles in (c). Lo- 
cal minima in conductance are found at the nearest-neighbor 
(NN) sites, (d) Evolution of tunneling spectra taken along 
the line indicated in (a) and (b) (corresponding to the a-axis 
or gap node direction), at different distances from position 
"A". Only one peak appears, corresponding to the impurity 
bound state at -3 mV. The peak energy does not depend on 
the measurement position. 



III. RESULTS AND DISCUSSION 

In most regions of the sample probed in this study, 
we observed isolated Zn resonance sites with the charac- 
teristic " cross shaped pattern^^ , as shown in Fig. [1] the 
concentration of these sites was approximately 0.45% per 
Cu atom (7^8 counts in a 15 x 15 nm^ field of view), 
which is comparable to the concentration of Zn in the 
bulk. 

According to theoretical worki^"— , when the distance 
between two Zn atoms becomes shorter than IOoq {oq is 
Cu-O-Cu distance), a spatial modulation of the LDOS 
near Ep should appear in the vicinity of the Zn atoms 
as a result of interference between the bound states. In- 
deed, we found some regions of the sample where such 
a LDOS modulation partially appears, as shown by the 
G{r, —2 mV) map for a 15 x 15 nm^ field of view in 



Fig. [U^a). Several bright islands are apparent, in which 
a two-dimensional (2D) LDOS modulation along the a- 
and 6-axes is visible with a period of approximately 5.4 
A (a/2 X ao); this is shown more clearly in Figs. [D^b) 
and (c). Figures[2l^e)-(g) show series of tunneling spectra 
measured at intervals along the three lines indicated in 
Fig. [2ta), which all pass through the bright islands. All 
spectra taken in the bright islands exhibit a NZEP corre- 
sponding to the impurity bound state. This implies that 
the 2D LDOS modulation is strongly related to the impu- 
rity bound state. However, the observed spatial pattern 
is a lattice- shaped pattern, which is significantly different 
to the cross-shaped pattern observed in the vicinity of an 
" isolated single^'' impurity as shown in Fig. [Ijb) . There- 
fore, it is conceivable that two or more impurities are 
present in the bright islands in the G{r, — 2mV) map, al- 
though the exact positions of these impurities could not 
be identified. If we assume that there are at least two 
impurity atoms in each bright island in Fig. [^Ja), there 
are more than 12 impurity atoms in this field of view, 
corresponding to a concentration of at least 0.8% per Cu 
atom. 

In addition to these lattice- shaped pattern, we observed 
several more features that are significantly different to 
those present near an " isolated single" impurity, as illus- 
trated by the spectral surveys on the 2D LDOS modu- 
lation in Figs, me)-(g). First, the energy of the NZEP 
exhibits a pronounced spatial variation. For example, in 
Fig. [UJe) the peak energy drastically shifts from -1 meV 
at site A to -4 meV at site B, despite the relatively short 
distance {'-^ 7 A) between these sites. Second, positive 
energy states appear, as shown in Fig. ^i) where the 
peak energy abruptly shifts from -2 meV to -1-5 meV at 
site C. Third, a splitting of the impurity bound state is 
observed. For example, the spectrum at the bottom of 
Fig. [2jg) has only one broad peak at positive energy. As 
one moves along Line 3 to site D, a negative energy peak 
gradually appears at -4 meV. At site D, two peaks at 
-4 meV and -1-3 meV are observed, and are shown more 
clearly in Fig. [2l^d). 

These three peculiar phenomena in the spectral sur- 
veys can be explained as a result of quantum interfer- 
ence among the impurity bound states by comparison 
with theoretical calculation o^^i^^'^^" — . The simplest sit- 
uation for which this interference effect can theoreti- 
cally be discussed is the two-impurity case. It is pre- 
dicted that at most four bound states will be generated 
(two states at negative energy and two states at positive 
energy) ^^'^"i^^" —. According to these calculations, the 
spectral weights and the energies of these states depend 
on the measurement position and the distance between 
the impurities, respectively. Furthermore, it has been 
predicted that (i) an abrupt change in the bound state 
energy, (ii) an appearance of positive energy states, and 
(iii) a split of the bound state can occur as a result of 
quantum interference when the impurity atoms are found 
in a particular configuration; these are all phenomena 
observed in our current study. Although a quantitative 
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FIG. 2. (Color) (a) Conductance G(r, Vb) map (field of view 15 x 15 nm'^) measured at a bias voltage of Vb = -2 mV. 
The white arrows indicate the crystal axes determined from the simultaneously obtained topographic image (inset), which was 
measured with Vb = 150 mV and a tunneling current of /sot ~ 600 pA. (b),(c) Magnified conductance map and topographic 
image for the region outlined by the blue box in (a). The positions of Bi atoms are indicated by black dots in (b) and red 
dots in (c). (d) Representative tunneling spectra measured at positions "A" to "D" indicated in (a). For clarity, the spectra 
at positions "A", "B" and "C" have been shifted upward by 0.75, 0.5, and 0.25, respectively. Blue and red triangles indicate 
peak positions at negative and positive energy, respectively, (e)-(g) Tunneling spectra measured at intervals along Lines 1-3 
indicated in (a). The four spectra in (d) are shown in white in (e)-(g). 



comparison of our experimental data with the theoret- 
ical calculations is not possible due to the uncertainty 
in the impurity positions, the observed features quali- 
tatively agree with the calculations. Because the three 
features ascribed to the quantum interference effect occur 
on the 2D LDOS modulation, we expect that the modu- 
lation also originates from the quantum interference. 

In order to confirm the origin of the 2D LDOS modula- 
tion, we investigated the energy evolution of the G(r, Vb) 
map and performed Fourier analyses as shown in Fig. |3l 
In the Fourier transform (FT) image of the G(r, — 2mV) 
map shown in Fig. EJ^g), the four broad spots indicated 
by red dashed circles correspond to the 2D LDOS mod- 
ulation. In the cut taken along the blue line in this FT 
image, a broad peak appears at g ~ (l/-\/2)(27r/ao), as 
indicated in Fig. [3jj). This is consistent with the real- 
space analysis regarding the modulation. As the energy 
is increased up to |£:| = elVel = 10 meV, the 2D LDOS 
modulation becomes less well defined and is no longer 
visible at energies above 14 meV. Instead, the conduc- 
tance map then shows an inhomogeneous pattern reflect- 
ing the inhomogeneous superconducting gap distribution, 
which commonly exists in cuprate superconductors^^—. 
The period of the 2D LDOS modulation does not change 
with increasing energy within our g-spacc resolution (^ 
0.025 X 27r/ao), as shown in Fig. 3(j). In addition to this 
non-dispersive nature, the amplitude of the LDOS modu- 
lation is strongly suppressed by taking the ratio between 
the conductance at a positive and a negative bias volt- 
age [i?(r, Vb) = G(r, -hVB)/G(r, -Vb)], as shown in Fig. 
[3]Jf) and the corresponding FT image in Fig. EJi). 



From these results, several possible explanations for 
the origin of the observed LDOS modulation can be ex- 
cluded. First, the non-energy-dispersive nature of the 
LDOS modulation and the suppression of amplitude in 
the i?(r,VB) map are inconsistent with the properties 
of quasiparticle interference (QPI), which has been ob- 
served for several cuprate superconductors^^"— , since the 
amplitude of the QPI pattern should be enhanced in the 
i?(r, Vb) map^^'^'^. Furthermore, the observed 2D LDOS 
modulation is completely different to the well-known non- 
energy-dispersive LDOS modulation with a period of ^ 
4ao X 4ao (so-called " checkerboard modulation''^ ) and ^ 
4/3aox4/3ao (so-called " nano- stripe") modulation that 
are observed in both underdoped cuprate o^^'"^^" — and 
around vortex cores^i^. In our case the observed modu- 
lation is aligned along the a- and 5-axes, which are tilted 
by 45° from the direction of the checkerboard and nano- 
stripe modulations. 



The features that we observe agree rather well with 
the predictions of calculations regarding the quantum 
interference effect of randomly distributed impurities^. 
These calculations indicate that 2D LDOS modulation 
can partially appear along the a- and 6-axes with a pe- 
riod of 5.4 A as a result of the interference of impurity 
bound states when the impurity concentration is 0.5% 
per Cu atom, which is almost the same as that in our 
sample. Thus, it is plausible that the 2D LDOS modu- 
lation and the three peculiar phenomena observed in the 
spectral surveys reflect the quantum interference effect. 
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FIG. 3. (Color online) (a)-(e) Conductance G(r, Vb) maps at Vb ~ -2, -6, -10, -14, and -30 mV for the same field of view as 
in Fig. [IJa). (f) Conductance ratio map i?(r, Vb) = G(r, -|-Vb)/ G(r, — Vb) at Vb = 2 mV. (g)-(i) Fourier transform images 
corresponding to maps (a), (e), and (f), respectively. The four black circles indicate the Fourier spots of the atomic lattice. 
In (g), four broad spots corresponding to the 2D modulation are indicated by red dashed circles; no such spots are visible in 
(h) and (i). (j) FT intensity taken along the cut indicated by the blue arrow in (g) for various values of Vb- The black line 
represents the corresponding line profile in the ratio map at Vb = 2 mV. 



IV. SUMMARY 



In summary, we have investigated the effect of mul- 
tiple impurities on the LDOS in Zn-Bi2212 using scan- 
ning tunneling spectroscopy. Our results reveal several 
peculiar features: (i) the 2D LDOS modulation vi/ith a 
period of approximately 5.4 A along a- and 6-axis; (ii) 
the abrupt change in the energy of the impurity bound 
states with position; (iii) the existence of positive energy 
states; (iv) the splitting of the impurity bound state. All 
of these features have been predicted in several theoreti- 



cal studies focused on the effect of quantum interference 
on the impurity bound states. Due to the good quali- 
tative agreement between our results and these theoret- 
ical calculations, we conclude that the observed features 
are fingerprints of the quantum interference of impurity 
bound states. Furthermore, one may expect that the ac- 
cumulation of these states at several energies induced by 
the interference creates the impurity band that is inti- 
mately related to the bulk transport properties. Thus, 
the quantum interference of impurity bound states ob- 
served in this study plays a crucial role in understanding 
how the microscopic impurity effect is linked to macro- 
scopic properties such as the suppression of Tc. 
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